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Abstract 



A new method is introduced for studying boundary value problems for a class of linear PDEs 
with variable coefficients. This method is based on ideas recently introduced by the author for the 
study of boundary value problems for PDEs with constant coefficients. As illustrative examples the 
following boundary value problems are solved: (a) A Dirichlet and a Neumann problem on the half 
line for the time-dependent Schrodinger equation with a space dependent potential, (b) A Poincare 
problem on the quarter plane for a variable coefficient generalisation of the Laplace equation. 



1 Introduction 



The main aim of this paper is to introduce a methodology for solving boundary value problems for 
a class of linear PDEs with variable coefficients. As illustrative examples the following boundary 
value problems will be solved: 

(i) A Dirichlet problem for the time-dependent Schrodinger equation on the half line, 

iQt + Qxx + u{x)q = 0, < x < oo, t > Q, (1.1) 

q{x,0) = qo{x), 0<x<oo, (1.2a) 

q{0,t)=go{t), t>0, (1.26) 

where q{x, t) is a complex-valued function, u{x) is a given real-valued decaying function, and qo{x), 
go (t) are given complex- valued functions with sufficient smoothness and decay which are compatible 
at x = t = 0, i.e. go(0) = 5o(0). 

(ii) A Neumann problem for the time-dependent Schrodinger equation, i.e. the problem defined 
by (1.1), (1.2a), and by the equation 

qx{0,t) = gi{t), t>0, (1.26)' 

where qo{0) = 51(0). 

(iii) A Poincare problem for a variable coefficient generalisation of the Laplace equation on the 
quarter plane, 

Qxx + Quy + u{x)q = 0, 0<x<oo, 0<y<oo, (1.3) 
gj^(x,0) +7ig(x,0) = /(x), < .x < 00, (1.4a) 

9a; (0,y)+/3gj/(0,2/)+ 725(0, y) = 5(y), 0<y<oo, (1.46) 

where q{x,y) is a rcal-valucd function, u{x), f{x), g{y), are given real-valued functions with ap- 
propriate smoothness and decay, and /?,7i,72 are given real constants. 



A. The Case of Constant Coefficients 

A new method for studying boundary value problems for linear and for integrable nonlinear 
PDEs in two independent variables was introduced in [1]. For linear PDEs in a convex polygon it 
involves two novel steps [2]: 

(a) Construct an integral representation of q in the complex k-plane. For example, for a linear 
dispersive evolution equation on the half line this representation is [3] 

1 1 r 

q[x,t) = — dke'^'^-'^'^^^^Qik) + — dke'^''-'''^^^^g{k,t), (1.5) 
27r 7_oo 27r Jqjj^ 

where ui{k) is the associated dispersion relation, qo{k) is the Fourier transform of (/(x,0), g{k,t) is 
an appropriate time transform of q{0, t) and of {dlq{0, t)}i~^ {N is the order of the highest spatial 
derivative), and dD^ is the boundary of oriented so that is on the left of dD^, where -0+ 
is the domain in the upper half of the complex /c-plane defined by Im uj(k) > 0. 

In the case of equation (1.1) with u{x) = 0, it can be shown that Lo{k) = ik"^ , dD^ is the 
boundary of the first quadrant of the complex A;-plane, and 

/•oo 

qo{k)= dxe-'''''qo{x), Im /c < 0, (1.6) 
^0 
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g{k,t) = kgo{k,t)-igi{k,t), k e C, 



(1.7) 



where go and gi axe defined by 




dTe"'"^qoiO,T), gi{k,t)= f 



(1.8) 



We note that in the Dirichlet problem g{k,t) involves the unknown function gi(k,t), while in the 
Neumann problem g{k,t) involves the unknown function go{k,t). 

In the case of the Laplace equation on the quarter plane, q{x, y) admits an integral representation 
defined along the oriented boundary of the first quadrant of the complex fc-planc [2], [4]. For the 
Dirichlet problem this representation involves the unknown Fourier transforms of qy{Q,x) and of 
ga;(0,y). 

(b) Use the invariant properties of a certain algebraic global relation to express the unknown 
transforms in terms of transforms of the given initial and boundary conditions. 

In the case of equations (1.1), (1-2) with u{x) = 0, it can be shown that this yields 



where q{k,t) denotes the Fourier transform of q{x,t). We note that the rhs of equation (1.9), in 
addition to the known functions c/q and qo, it also involves the unknown function q{—k, t). However, 
this term does not contribute to q{x,t); indeed it gives rise to the term 



which vanishes according to Jordan's lemma, since both exp[zA;a;] and q{—k,t) axe analytic and 
bounded in the first quadrant of the complex /c-plane. Thus dropping the last term of the rhs of 
equation (1.9) and substituting —igi in the rhs of equation (1.7) we find 



which involves only transforms of the known functions qo{x) and go(x). 

In the case of equations (1.3), (1.4) with u{x) = 0, the analysis of the global relation implies that 
the unknown transforms can be computed through the solution of a scalar Riemann-Hilbert problem 
[4]. For particular values of /3, 71, 72 (which include both the Dirichlet and the Neumann cases) this 
Riemann-Hilbert problem can be bypassed and the unknown transforms can be computed using 
only algebraic manipulations (just like the case for equations (1.1), (1.2)). 

For evolution PDEs there exist at least three different ways of implementing step (a). Use: (i) 
The spectral analysis of the associated Lax pair [2]. (ii) A reformulation of Green's theorem [5]. 
(iii) The deformation of the Fourier transform representation [6]. Methods (i) and (ii) can also be 
used for elliptic PDEs. 

For particular PDEs, such as equation (1.1) with u{x) = 0, it is possible to avoid the construc- 
tion of an integral representation in the complex A;-plane, and to use instead the usual integral 
representations along the Re /c-axis. However, in general the complex fc-plane cannot be avoided; 
in particulax this is the case for an evolution PDE involving a third oxder dexivative. 

B. The Case of Variable CoefRcients 

The analysis of certain classes of PDEs with variable coefficients is also based on the two novel 
steps (a) and (b) mentioned in A above. 



-igi{k,t) = kgo{k,t) - qo{-k) + e''' ^qi-k,t), lmk>0, 



(1.9) 




g{k,t)=2kgo{k,t)-qo{-k), lmk>0, 
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(a) The construction of an integral representation. For both evolution and elliptic PDEs this can 
be achieved using the simultaneous spectral analysis of the associated Lax pair. A Lax pair for 
equation (1.1) is 

fJ'xx + (uix) + fc^)/x = q (1.10a) 

fit + ik'^ H = iq. (1.106) 

Indeed, equations (1.10) are compatible iff q{x,t) satisfies (1.1): Applying the operator dt + ik'^ to 
equation (1.10a), noting that this operator commutes with 5^ + u{x) + fc^, and using (1.10b) we 
find (1.1). Similarly a Lax pair for equation (1.3) consists of equation (1.10a) together with 

Hyy - k^n = -q. (1.11) 

For evolution PDEs, an alternative approach to constructing integral representations is to use a 
completeness relation of the associated space dependent eigenf unctions. This is the analogue of the 
deformation of the Fourier transform, see (iii) in A. For equation (1.1) the relevant completeness re- 
lation is the restriction to the half line of the classical relation associated with the time-independent 
Schrodinger equation [7]: Let u{x) G -Lf, i.e. assume that dx{l + x'^)\u{x)\ exists. Let ijj{x,k) 
and (f){x, k) be the following solutions of the time- independent Schrodinger equation 

Xxx + {u{x) + k')x = Q, x>0, fceM: (1.12) 

1 1"°° 

tp{x,k) = e''''' + - dS,smk{x - OuiO''Pi^,k), x > 0, lmk>0, (1.13) 

(f)(x, k) = e-^'^^ --[ d^ sin k(x - OuiOH^, k), x >0, k e C. (1.14) 

k Jo 

Let kj = ipj, j = 1, ■ ■ ■ , n,pj > 0, be the zeros of a(k), where a{k) is defined by 

a{k) = ^{Hx-i^^x), ImA;>0. (1.15) 

Then 

where d{k) denotes the derivative of a{k). 

(b) The global relation. This equation can be obtained using the formal adjoint of the given PDE. 
For example, if (1.1) is valid in a simply connected domain with boundary F, the global relation is 

J^{{qQ)dx 

X 

)dt] = 0, (1.17) 

where Q is any solution of the "adjoint" equation 

-iQt + Qxx + uQ = 0. (1.18) 
Indeed, if q and Q satisfy equations (1.1) and (1.18) then 

{qQ)t + i{qQx - qxQ)x = o, 

and Green's theorem yilcds (1.17). 

For equation (1.13) the global relation is 

J^{M - qQx)dy + {qQy - qyQ)dx'^ = 0, (1.19) 

where Q is any solution of equation (1.3). 
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C. Organization of the Paper and Main Results 

Using the completeness relation (1.16) as well as the global relation (1.17), it will be shown 
in section 2 that the solution of the initial-boundary value problem (1.1), (1.2) is given by the 
following formulas: 



27r 



V'(x, kj 



-ik-'t 



a{k) 



Qo{k) 



q{k,t) -iY] 
qo{k) 2ka{k) 



+ 



a{k. 
9o{k,t) )■ , 



(1.20) 



(1.21) 



^'(O,^) V(0,A;)- 

where il^ix, k), (f){x, k), a{k) are defined in terms of u{x) by equations (1.13)-(1.15), goik, t) is defined 
in terms of the boundary condition q{0,t) = go{t) by equation (1.8a), and qo{k), qo{k) are defined 
in terms of the initial condition qo{x) by 

poo poo 

Qoik) = / dx(l){x,k)qo{x), qoik) = / dxil){x,k)qo{x), kER. (1-22) 
Jo Jo 

The solution of the Neumann problem is given by similar formulas but q{k,t) is now given by 



q{k,t) 



-ik'^t 



qo{k) + ik 



Qoik) 



+ 



2a{k) 



MO,k) MO,k) 



■9i{k,t) > , 



(1.23) 



where gi is defined in terms of the boundary condition qx{0,t) = gi{t) by equation (1.8b). 

It is possible to deform some of the terms appearing in the integral representation of q{x, t) 
to an integral along the boundary of the first quadrant of the complex /c-plane. For example it is 
shown in section 2 that equations (1.20), (1.21) imply 



+- 



1 



%ik) 
' a{k) 



+ 2kgo{k,t) 



(1.24) 



27r Jqd, ^(0, k) 

where dDi denotes the boundary of the first quadrant of the complex A;-plane encircling the points 
{kj}i, from the right, see Figure 1.1. 

Im k 



kn*) 



kii) 



Re k 

Figure 1.1: The contour dDi associated with equations (1.1), (1.2). 
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If u{x) = 0, then tp = exp{ikx), (p = ex.p[—ikx), a{k) = 1, and the representation (1.24) coincides 
with the representation (1.5). The representation (1-24) is convenient for analyzing the long time 
behavior of q(x,t), see section 2. 
The particular example of 



2^ 

cosh^ p{x — Xq) 



^(^) = 7^-7 7' P>^^ XQ>0, (1.25) 



is also discussed in section 2: in this case ■ijj,(l),a, can be computed cxphcitly. 

In section 3, using the spectral analysis of the Lax pair (1.10a), (1.11), as well as the global 
relation (1.19), we will solve equations (1.3), (1.4), see propositions 3.2, 3.3. 

In section 4 we discuss further the above results. In the appendix, for the sake of completeness, 
we present a simple derivation of equation (1.16) (see also [7]-[9]). 

2 The Time-Dependent Schrodinger Equation 

(a) An Integral Representation 

The completeness relation (1.16) yields 

where 

poo 

q{k,t)= dxq{x,t)(f){x,k), k eR. (2.2) 
Jo 

If q(x,t) satisfies equation (1.1) and (p{x,k) is defined by (1.14) then 

(e'^^cp)^ + ^e^''*(#x - qx<f>)x = 0. (2.3) 
Multiplying equation (2.2) by e**^^*, differentiating with respect to t, and using equation (2.3) we 
{e^'^'qik^t))^ = ie"'"\q{0,t)MO,k) - q,{0,t)(t>{0,k)). 



find 
Thus 



e"'"'q{k,t)=q{k,0)+iMf^,k)go{k,t)-i<f>{0,k)gi{k,t), (2.4) 

where go and gi are the t-transforms of q{0,t) and of qx{0,t) defined by equations (1.8). 
In summary, q{x,t) is given by equation (2.1) where q{k,t) is defined by equation (2.4). 

(b) The Global Relation 

Let ijj{x,t) be defined by equation (1.13), and let 

/•oo 

q{k,t)= dxq{x,t)ip{x,k), ImA;>0. (2.5) 
Jo 

In analogy with equation (2.4) we find 

e'''^^q{k,t) = q{k,0) + iilJa:{0,k)go{k,t)-iip{0,k)gi{k,t), lmfc>0. (2.6) 



6 



1. The Dirichlet Case 

Solving equation (2.6) for gi and substituting the resulting expression in (2.4) we find 



It is important to note that the term q{k,t) and the associated terms q{kj,t) do not contribute to 
q{x,t). Indeed these terms give rise to the terms 

The function ■^(a;, k)e~'^^^ is analytic for Im A; > 0, thus since x > 0, the function ip{x, k) is also 
analytic for Im/c > 0. Hence all the terms appearing in the integrant of the first term in (2.8) are 
analytic for Im /c > 0. Furthermore, ^(0, A;) 7^ for Im A; > [8]. Also the contribution from the 
zeros of a{k) is cancelled by the second term of equation (2.8). Hence Jordan's lemma applied in 
the upper half of the complex A;-plane implies that the term (2.8) vanishes. 

In summary, q{x,t) is given by (2.1) where q{k,t) is given by equation (2.7) without the first 
term of the rhs of (2.7). 



2. The Neumann Case Solving equation (2.6) for go and proceeding as above we find 



-ikh f ^ /,x MO,k) 2ka{k) o 



Since (p{x,k) satisfies (1.14), (p{0,k) = 1, (px{0,k) = —ik. 
2.1 The Completeness Relation on the Full Line 

There exist a particular class of functions u{x), —00 < x < 00, called reflectionless potentials. 

For the restriction of this class of potentials on the half line, it is more convenient to use the 
completeness relation associated with the full line instead of equation (1.16), see the Appendix. 
This relation is defined by equation (1.16) where (j){x, k), a{k), are replaced by ^{x, k), A(k), which 
are defined as follows: 

$(x,A;) = e-*'^^ - ^ / d^smk{x-Ou{0^^,k), keR, (2.10) 

A{k) = ^{^ijx-^xiP), lmk>0. (2.11) 

The functions $e**^^, A{k) arc still analytic for Imk > 0, furthermore tplx, k) is analytic for Im/c > 
provided that x > 0. Thus the arguments used for the derivation of g(x, t) apply mutatis-mutandis 
if ^, a are replaced by ^,A. In summary: 

Proposition 2.1 Let the real-valued function u{x), —00 < x < 00, satisfy dx(l-|-x^)|u(x)| < 
00. Define ■ip(x,k), ^{x,k), A{k), in terms of n(x) by equations (1.13), (2.10), (2.11). Let kj = ipj, 
Pj > 0, j = 1, • • • , n, be the zeros of A{k) for Im fc > 0. The solution of equations (1.1), (1-2) is 
given by the following formulas: 

= ± -,f:i|M,ft,,), (2.12) 
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where qq is defined in terms of the boundary condition q{0,t) = go{t) by equation (1.8a), while 
Qo(,k), qo{k) are defined in terms of the initial condition qo{x) by equation (1.22b) and by 



go 

For the Neumann case 



ik)= dx^{x,k)qo{x), keR. (2.14) 
Jo 



«M) = e-" {m) - 2^«(*) + |^S.(M)} . (2.15) 

Example 2.1 Let u{x) be given by equation (1.25). In this case A{k), ipix^k) are defined by 
equations (A. 13), <t{x,k) = a{k)ijj{x,k), and ki = ip. Thus 

^(^i) = 7r^> ki) = — ^, ^{x, ki) 



2ip^ ' 2coshp(x — xo) ' ' 2coshp(x — xq) 

Hence, equations (2.12), (2.13) simphfy to the following expressions: 



q{x, t) = ^(^J dkij{x, k)Qik, t) + ^i(x)Qi(t)^ , (2.16) 



Q{k, t) = e-'^ * <^ dxi^ix, k)qo{x) - dxi;{x, k)qo{x) + -^^^go{k, t) \ , (2.17) 



coshp(x — xq) 

2.2 Alternative Representations 

Inserting equation (1.21) into (1.20) we find 



^i(^) = u / Qlit) = e^'* / dxMx)(loix). (2.18) 



7^ 



-e J 



27r7_oo V'(0,A;) a{k) ^ ^{Q,kj) a{kj) \ 



The term 

ip{x,k) ikH(lo{k) 



-e 



V'(0,A;) a{k) 

is bounded and analytic for 7r/2 < argfc < tt. Thus using Cauchy's theorem the integral of this 
term along the negative real axis can be rotated to an integral along the positive imaginary axis; the 
associated residue sum can be incorporated by encircling the poles from the right. The integrant 
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of the last integral of the rhs of (2.19) is also bounded and analytic for 7r/2 < argA; < tt. Indeed 
iIj{x, k)/tp{0, k) is analytic for Irak > and 

Jo 

thus this term is bounded for k G [7r/2, tt] U [37r/2,27r]. Hence the integral of this term along the 
negative real axis can also be rotated to an integral along the positive imaginary axis. Hence 
equation (2.19) becomes equation (1.24). 



2.3 The Long t Asymptotics 

Suppose that equation (1.1) is valid for < t < T. Then go{k,t) in (1.24) can be replaced by 
go{k,T). Indeed, the difference of these terms gives the contribution 

which vanishes due to analyticity considerations. In particular, if goit) decays as t ^ oo, go{k,t) 
can be replaced by go{k), 

roo 

go{k) = / dre"" ^go{r). (2.20) 
Jo 

Equation (1.24) with go{k,t) replaced by go{k) has the advantage that it involves time-dependence 
only in the explicit exponential form exp(— zA;^*) and exp(— iA;|t). This formulation is most conve- 
nient for the analysis of q{x,t) as t ^ oo, see [10]. 



3 The Generalized Laplace Equation 

We first derive an appropriate integral representation for equation (1.3), as well as the relevant 
global relations. We then apply these formulae to the solution of the boundary value problem 
(1.3), (1.4). 

For the sake of simplicity we assume that the function a{k) has no zeros for Im A; > 0. 



Proposition 3.1 Assume that there exists a real-valued function q{x, y), < x < oo, < y < oo, 

which satisfies equation (1.3). Assume that q{x,y) has sufficient decay as a; ^ oo, y oo, and 
it also has sufficient smoothness all the way to the boundary. Then q{x,y) admits the following 
integral representation: 

\ roo 

y) = ~—J^ dke-^y{p{k)il){x, k) - p{k) 4j{x, k)) 

where the functions go{k), p{k), are defined in terms of the boundary values q{x,0), q{0,y), by 

roo 

goik)= dye''yq{0,y), Re A; < 0, (3.2) 

Jo 
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foik) 



dxq{x, 0)(p{x, k) - 
V'iOjA;) 



/•oo 

/ dxq{x,{))'4){x,k), ImA;>0. (3.4) 
Furthermore, the following global relations are valid 

/ dx[kq{x,Q) - qy{xM<i>{x,k) + d?/e^^[(z(0, y)(/)^(0, fc) - ^0.(0, y)0(O, /c)], Re A; < 0, (3.5) 
Jo Jo 

roo poo 

/ dx[kqix,0)-qyix,0)]i;ix,k)+ dye*^^[g(0, y)V'^(0, /c) - (/^(O, y)V'(0, fc)], Im A; > 0. (3.6) 
Jo Jo 

Proof 

We will construct four functions lJ:i{x, y,k), - ■ ■ , /X4(x, y, k) which satisfy both equations (1.10a), 
(1.11), and which for < x < oo, < y < oo, are bounded and analytic in k for k in the first, 
■•• , fourth quadrant of the complex fc-plane. Equations (1.10a), (1.11) are invariant under the 
transformation k —k, thus it is sufficient to construct only the functions ni and /X2, see Figure 
3.1. These functions are given by 



2kiii{x,y,k) = h{x,y,k) + ^^^'^^ /2(y, fe), < argfe < -, 

llJ f X ki I IT 

2kn2ix,y,k) = Ii{x,y,k) - ^i^^^—^l2{y,-k), - <argA;<7r, 



(3.7) 
(3.8) 



where 
Ii{x,y,k) 



ia{k) 



POD 

4>ix,k) / d^qi^,y)2l;{^,k) 

J X 



+ 



i^{x,k){ 1^ diq{i,y)<i>{i,k)-^^^ diq{i,y)^{i,k) 



, ImA;>0, (3.9) 



l2{y,k)= dr/e-'=(^-'')g(0,r/) + / dr]e^^y-i\{Q,r]), Re )fc > 0. (3.10) 

Jo Jy 

Indeed, a particular solution of (1.10a) which is analytic for lm.k > 0, is given by (see the Appendix) 



2kiJ, = — 
ia 



' / diq{i, y)V(C, fc) + V / diq{i, y)0(^, k) 
Jx Jo 



We add to this particular solution the homogeneous solution p{k,y)ip{x,k) and we choose p{k,y) 
so that at X = we have an identity: 



2kn{x, y, k) = h{x, y, k) + 
We now use equation (1.11) to determine ij,{0,y,k): the functions 

niO, y,k) = ^l2{y,k), n{Q,y,k) = -^hiy^-k), 



(3.11) 



(3.12) 



are particular solutions of equation (1.11) evaluated at x = 0, and are analytic for Re > 0, 
Re/c < 0, respectively. Substituting equations (3.12) into (3.11) we find (3.7), (3.8) 
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2kiJ,2{x,y,k) „ 2kfii{x,y,k) 



2kn'2.{x,y,-k) 2kfii{x,y,-k) 



Figure 3.1: The domains of analyticity of the functions /xi and ^2- 

We can use the functions appearing in Figure 3.1 to formulate a Riemann-Hilbcrt problem. 
To this end we need the large k behavior of /Hj, j = 1, 2, as well as the relevant "jumps". The 
definitions of Ii, I2, imply 



fe/xi {x,y,k) = O [ -j^ , k ^ 00 in 0< arg k < —; 
kiJL2 {x^y^k) = O [ -j\ , k ^ 00 in — < arg k < tt. 



(3.13) 



In order to compute the jump across the ray arg A; = 7r/2 we subtract equations (3.7), (3.8), and 
use the definition of I2, 



2kixi{x,y,k) - 2kn2{x,y,k) 



■(/'(x, k) 



e-''ygo{k)+e''ygo{-k) , arg A; 



TT 



(3.14) 



In order to compute the jump across the ray arg A; = we note that since both iJ,i{x,y,k) and 
IJ,2{x,y, —k) satisfy equations (1.10b), (1.11), it follows that their difference satisfies the homoge- 
neous versions of (1.10b), (1.11). Also exp[A;y] cannot appear in this equation since it is unbounded. 
Hence 

2ki^i{x,y,k) - 2kii2{x,y,-k) = e~'^y[pi{k)ij;{x,k) + p2{k)ijj{x,-k)], arg A; = 0; (3.15) 
letting y = 0, x ^ 00, we find 



pi{k) = p{k), p2{k) = -p{k). 

Let us denote the rhs of equations (3.14), (3.15), by Ji2(A;), JiA{k) respectively. The solution 
of the above Riemann-Hilbert problem is 



2kp{x,y,k) = ^ 



" ^'^^ 

I — k 



l-k 



or 



l+k l-k 



(3.16) 
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where fi = ni for < arg k < ^ and ^ = ^2 for 7r/2 < arg k < n. 

Equation (1.10a) implies q = limjt_»oo(/c^/x); this equation and equation (3.16) yield 



1 r°° If 



dlJuil), 



which is equation (3.1). 

If r is the boundary of the first quadrant of the complex z-plane, the global relation (1.19) 
becomes 

/ dx[q{x,0)Qyix,0)-qy{x,0)Qix,0)]+ / dy[q{0,y)Q^iO,y) - q^iO,y)Q{0,y)] = 0. 
Jo Jo 

Letting 

Q = e^y<t){x, k), Q = e^yi){x, k), 
in the above equation we find equations (3.4), (3.5). QED 

Remark 3.1 The rhs of equation (3.1) is a real-valued function: The integrant of the first integral 
is imaginary, also i/j{x,k) is real for k purely imaginary thus the integrant of the second integral is 

imaginary as well (due to the term dk). 

The integral representation (3.1) and the global relations (3.4), (3.5) can be used to solve a 
variety of boundary value problems. As an example we consider the boundary conditions (1.4). 

Proposition 3.2 The solution of the boundary value problem (1.3), (1.4) with /3 ^ is given by 
equation (3.1) where the functions p{k) and goik) can be computed in terms of the given functions 
f{x) and g{y) as follows: The sectionally analytic function {go(,k), go{—k)} satisfies the scalar 
Riemann— Hilbert problem, 

go (k) is analytic for Re A; < 0, 
go{k) = o(^j, k^oo, 



go{k) + mgo{-k) = ^^tm_q^o,0)+F{k), argA; = f , 

[ go(-fc) + J{-k)goik) = ^,f,%% qiO, 0) + Fi-k), argA. = f, 



(3.17) 



where 



''^^^ = fc-7i^.(0,fe) + (72-W(0,fc) ' I'(k)=MO,k) + ij,-pk)mk), (3.18) 



p. , ^ _ (fc-7i)g(fc) + (fc + 7i)iJ(-fc) 

ik-^i)D{k) ' ^^-'^^ 

roo roo 

H{k)= dxfix)i;{x,k)+i^iO,k) dye^ygiy). (3.20) 
Jo Jo 

The function p[k) is given by 

"W-I^^ + ^W. 0<a.g*<|. (3.21C) 

^'(O, k) dxf{x)(t){x, k) - (pjO, k) Jq°° dx/(x)V'(x, k) /qoim 
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Proof 

Solving the boundary conditions (1.4) for qy{x^ 0), qx{^, y), substituting the resulting expressions 
in (3.5), and integrating by parts the term involving qy{0,y), we find 

roo 

{k + ji) dxq{x,0)(t>{x,k) + d{k)go{k) = P(t){0,k)q{0,0) +h{k), Re A; < 0, (3.22) 
Jo 

where the known functions d{k), h{k) are defined by 

roo roo 

d{k) = <t)x{Q, k) + (72 - /3fe)0(O, A;), h{k) = / dxf{x)<t>{x, k) + </.(0, k) / dye'^ygiy). (3.23) 

Jo Jo 

Similarly, equation (3.6) implies 

(k + 'yi) dxq{x,0)i;{x,k)+D{k)go{k) = l3ip{0,k)q{0,0)+H{k), -<argfc<7r, (3.24) 
Jo 2 

where the known functions D{k), H(k) are defined by (3.18b), (3.20). Both equations (3.22), (3.24) 
are valid for 7r/2 < argA; < tt. Solving these equations for /o(A;) we find 

f (u\ - ^ika{k)go{k) ^ ilj{0,k) dxf{x)i){x,k) - (l){0,k) dxf{x)ip{x,k) vr 

^'^"^ - mrnkT^) ^ mWT^) , 2 <arg/c<7r. 

Taking the complex conjugate of this equation and letting k —k we find an expression for fo{k) 
valid for < argA; < 7r/2. Substituting this expression into equation (3.3) we find (3.21a). 

We now show that equation (3.24) implies a Riemann-Hilbert problem for {go{k), go{—k)}. 
Taking the complex conjugate of equation (3.24) and letting k — > —k, we find an equation valid for 

< argA: < 7r/2, 

(-A; + 7i) / dxq{x,0)tp{x,k) + D{-k)goi-k) = Pi;{0,k)q{0,0) + H{-k),0<argk<Tr/2. (3.25) 
Jo 

Eliminating from equations (3.24), (3.25) the integral involving q(x,0)tf:{x,k) we find equation 
(3.17a). This equation defines the relevant "jump" across the positive imaginary axis; letting 
k —k we find the jump across the negative imaginary axis, i.e. equation (3.17b). The definition 
of go{k) implies go{k) = 0{l/k), A; — oo. QED 

Remark 3.2 It was shown in [4] that for certain particular cases of Poincare boundary conditions 
it is possible to bypass the relevant Riemann-Hilbert problem. Such cases were called algebraic. 
In what follows we discuss an algebraic case for equation (1.3). 

Proposition 3.3 The solution of the boundary value problem (1.3), (1.4) with /3 = 0, 71 < 0, is 
given by 



1 f /"O ri<x> 

q{x,y) = { dke-^yi){x,k)N{k)+ dke-^^i^ix^k)- 

2«7r yj^ Jo 1 



F{k) 



where N{k) and F{k) are defined in terms of the given functions f{x) and g{y) by equations (3.21b) 
and (3.19). 
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Proof The rhs of equation (3.1) can be rewritten in the form 



9o{k) 



dke-^yipix,k)p{k)+ / dke-^y^x^k) , , 

Jo V'(0,fc) 

+ — / dke''y^{x,k)p{-k)+ / dke^^i^ix^k) 



go(-fc) 



For P = 0, equation (3.17a) becomes 



9o{k) + ^^^go{-k) = F{k), argk = ^. 



(3.27) 



(3.28) 



In equation (3.27) we replace p{k) by equation (3.21a); also we express go{k) appearing in the first 
bracket of equation (3.27) in terms of go{—k) (using (3.28)), and go{—k) appearing in the second 
bracket of equation (3.27) in terms of go{k) (using again (3.28)). This yields the rhs of equation 
(3.26), and the additional terms 



goj-k) fc + 71 
V'(0, A;) - 71 



— / dke^y^x, 



k) 



gojk) fc-7i 
^'(O, A;) + 71 ' 



(3.29) 



where dDi, dD2 denote the boundaries of the first and second quadrants of the complex fc-plane. 
These terms vanish due to the analyticity of the relevant integrants. QED 



Remark 3.3 If 71 > 0, there exists an additional contribution due to the pole in the expressions 
(3.29). This contribution can be computed explicitly by evaluating the global relations at fc = —71. 
The procedure is similar to the one discussed in [4] . 



4 Discussion and Generalisations 

In this paper we have introduced a new method for solving boundary value problems for a class 
of linear PDEs with variable coefficients, in two independent variables. This method which can 
be used to analyze both evolution PDEs (such as equation (1.1)) as well as elliptic PDEs (such as 
equation (1.3)), is based on two novel steps: (a) Construct an integral representation for the solution 
q in the complex fc-plane, see for example equations (1.24) and (3.1). This representation involves 
appropriate transforms of boundary values of q (and of its derivatives); some of these boundary 
values are not prescribed as boundary conditions, (b) Use certain algebraic global relations to 
determine the transforms of the unknown boundary values; for evolution equations this step involves 
only algebra, while for elliptic PDEs, depending on the nature of the boundary conditions, it involves 
either algebra (see for example proposition 3.3), or the solution of a Riemann— Hilbert problem (see 
for example proposition 3.2). 

The remarkable fact about the new method is that the analysis of the global relations for the 
case of PDEs with variable coefficients, is similar to the analysis of PDEs with constant coefficients. 
This implies that such PDEs can be solved with the same level of efficiency as the corresponding 
PDEs of constant coefficients, except that one now uses appropriate base functions (for example 
■^(x, fc), 0(x, fc)) instead of the exponential function. For particular cases of variable coefficients 
these base functions can be constructed explicitly. 

For both evolution and elliptic PDEs the relevant integral representations can be constructed by 
performing the simultaneous spectral analysis of the associated Lax pair (see for example Propo- 
sition 3.1). This provides a generalisation to PDEs with variable coefficients of the formalism 
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introduced in [1]. For evolution PDEs tlicrc exist also an alternative procedure for deriving integral 
representations; this is based on the associated completeness relation and on contour deforma- 
tion (see for example Section 2). This approach generalizes to the case of variable coefficients the 
formalism introduced in [6]. 

In this paper for economy of presentation we have assumed that the solution of the given 
boundary value problem exists. However, it is possible to present rigorous theorems without the 
apriori assumption of existence; for the case of constant coefficients this is discussed in [3] and [11]. 

We conclude with some remarks. 

1. Both the Dirichlet and the Neumann problems of equation (1.1) can be solved using suitable 
completeness relations instead of the completeness relation (1.16). For the Dirichlet problem the 
relevant relation is 

^ > 2ttJ_^ a{k) V mk)J Y "(^^■) V V'(0,%); ^ ' 

This relation can be derived by performing the spectral analysis of equation (1.1a) with the condition 
/i(0, k) =0 [9]. Using (4.1) it follows that instead of equation (2.4) we now find 

e"''*q{k,t)=q{k,0) + kgo{k,t). 

Thus we find q{k,t) immediately without the need to use the global relation. Equation (4.1) is 
the generalisation of the sine-transform, while equation (1.16) is the generalisation of the Fourier 
transform (restricted on the half- line). It has been emphasized by the author thai lh(;re does not 
exist an analogue of the sine transform for problems with derivatives of third order. For example 
the spectral analysis of the equation 

^ + Xf = 0^ < X < oo, with /(O, A) = 0, 

in Rietz is incomplete [12]. An advantage of the new method introduced by the author is that it 
does not rely on the spectral analysis of the associated differential operator obtained by separation 
of variables. Similar considerations apply to PDEs with variable coefficients. Consider for example 
the Dirichlet problem for the PDE 

qt + qxxx + u{x)q = 0, < X < oo. (4.2) 

For this equation there does not exist an analogue of equation (4.1), but using the spectral theory 
developed in [13], [14], it is possible to obtain the analogue of equation (1.16). This implies that 
equation (4.2) can be solved using the method introduced in this paper, but cannot be solved by 
performing the spectral analysis of 

fJ-xxx + (u{x) + k^)fJ- = q{x), < X <oo, 

and demanding //(O, k) = 0. 

2. We note that the classical method of separation of variables fails for the boundary value 
problem (1.3), (1.4). Indeed, although both the PDE and the domain are separable, the boundary 
conditions at .x = are not separable. For example, in the particular case of 71 = the boundary 
condition at y = is consistent with the cosine transform, 

roo 

q{x, k)= I dyq{x, y) cos ky. 
Jo 
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The cosine transform of equation (1.3) yields 

qxx + u{x)q - k^q = f{x). 
But the cosine transform of the boundary condition at x = leads to the non-separable equation 

/•OO /"CXD 

^a;(0,fc) + 72g(0,fc) +/3 / dyqy{0, y) COS ky = dyg{y) cosky. 

Jo Jo 

3. Using the method introduced in this paper it is possible to analyze a large class of other 
boundary value problems. As an example we mention the following generalisation of the Helmholtz 
equation on the quarter plane, 

Qxx + Qyy + (xq + iu{x) + v{y))q = 0, < x < oo, 0<y<oo, (4.3) 

where a is a real constant and u{x),v{y), are given real- valued functions with sufficient decay as 

X — > OO, y ^ OO. 

4. The implementation of the new method to linear evolution PDEs with constant coefficients 
in two space dimensions is carried out in [6] . In this case one constructs an integral representation 
of q{xi,X2,t) in the complex {ki, A;2)-planes instead of the complex A;-plane. Using the results of [6] 
together with the methodology introduced in §2, it should be possible to solve boundary value prob- 
lems for certain classes of evolution PDEs in two space dimensions with variable coefficients. The 
main difference is that, in order to derive an integral representation for q{xi,X2,t), one must now 
use an appropriate completeness relation instead of the two-dimensional Fourier transform. In re- 
cent years progress has been made for constructing such completeness relations for two dimensional 
eigenvalue equations, see for example [15], [16]. 

Appendix — Completeness Relations 

Given u{x), define tp{x, k), (f){x, k), a{k), by equations (1.13)-(1.15). Evaluating the rhs of equation 
(1.15) at X = 0, and as a; — ^ oo, it follows that 

1 f'°° 1 / \ 

a(fc) = l + — j^ d^e^'^^ui^cPi^, k) = ^ (Vx(0, k) + ik^ k)) ■ (A.l) 

The definition of ijj{x, k) implies 

V'(x, k)e-'^'' = '^^2ikJ ^^V-~ e2^'=(«-^)j u{0^{^, k)e-'^^ , (A.2) 

thus V'(2^) k)e^''^^ is analytic for Im k > 0. The function (j){x, k) is an entire function of k; writing 
an equation analogous to equation (A.2) for the function 0(a;, k)e^^^ it follows that this function is 
bounded for Im/c > 0. Equations (A.l) imply that a{k) is also analytic for Imfc > 0. In summary 

V'(x, fc)e-^'=^ (l){x,k)e'^\ a{k) {A.Z) 

are bounded and analytic for Im > 0. 
We now consider equation 

IJ'xx + {u{x) + k'^)lJi = q{x), < x < OO, /c G M, (^.4) 
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where q{x) has sufficient smoothness as well as sufficient decay as x ^ oo. Using variation of 
parameters it follows that for A; G M, a particular solution of equation (A.4) is 

2ikaik)^ip = -(/)(x, k) d^qiOHC, k) + ^{x, k) ^ d^q{Cm, k). 

We choose the limits of integration in such a way that the rhs of the above equation is analytic for 
Imfc > 0: 

poo px 

2ika{k)iJ,+ {x,k) = ct>{x,k) d^qiOi'i^,k)+i;ix,k) d^qiOH^,k), lmfc>0. (A5) 

Jx Jo 

Indeed, the first term of the rhs of (A.5) equals 

/•oo 

{^{x,Oe"'n / d^e'''(^--\{Om,k)e-''^), 

Jx 

and exp[7;A;(^ — x)] as well as the two terms appearing in the parentheses are analytic for ImA; > 0, 
see (A. 3); similarly for the second term of the rhs of equation (A.5). 

Both functions n~^{x,k) and //+(x,— /e) satisfy the same equation (A.4), thus their difference 
satisfies the homogeneous version of equation (A.4), hence 

2ikfj,'^{x,k) -2ikii'^{x,-k) = p{k)ip{x,k) + p{-k)ilj{x,-k), keR. {A.6a) 

Evaluating this equation as x — oo, it follows that 

P^^^ = ^ /"«e)'/'(e,A:)- (A65) 
a{k) Jo 

Equation (A.5) implies that /x+(x,A;) is a meromorphic function of k for Im A; > 0. Also using 

the definitions of cp and -0 it follows that 2ikfi~^{x,k) = 0{l/k) as k ^ oo. These facts together 
with equation (A. 6a) define a Riemann-Hilbert problem for the sectionally meromorphic function 
{2iA;^+(x, A;), 2ikfi^ (x, —A;)}. Before solving this problem we need to compute the relevant residues. 
Let kj = ipj, pj > 0, j = 1, - ■ ■ , n, be the zeros of a{k) for Im > 0; the residue of 2ikfi'^{x, k) at 



k — kj IS 



ct>{x, kj) / d^qiOm kj) + 0(e, kj) / d^qim^, k 
Jx Jo 



kj) /o°° d^qiOH^^ kj 



{A.7) 

where we have used (p{x,kj) = Cjilj{x,kj) (which follows from equation (1.15), since a{kj) = 0). 
Equation (A. 6a) can be rewritten as 



2^kp^x, k)-J2t^j- |2zfc/x+(x, -/^) - E + 



pik)Hx, k) + pi-k)4,ix, -k) - x: ^ + f: 

1 1 

where Aj{x) denotes the rhs of equation (A.7). This equation imphes 
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Using q = \imk-*oo{k^ n) , equation (A. 8) yields 

1 ^ ^ 

qix) = 27r y dkil){x,k)p{k) -i^Aj{x), 

which is equivalent to equation (1.16). 

Remark A.l (j){x,k),il){x,k) are related by 

(l){x,k) =a{k)ip{x,-k) +b{k)ip{x,k), keR, {A.9) 

where 

2ikb{k) =il;{x,-k)(f)x{x,k)-4>{x,k)iljx{x,-k), keR. (AlO) 
The functions a{k),b{k) satisfy 

a{k)a{-k) - b{k)b{-k) = 1. (All) 

Since u{x) is real, 

ipix, -k) = ij{x,k), a{-k) = a(I), b{-k) = b(fij. {A.12) 

Remark A. 2 It is sometimes convenient to consider extensions of u{x) from < x < oo to 

— oo < X < oo, see the example below. The completeness relation for x on the full line is also given 
by equation (1.16), but with 0(x, k), a{k) replaced by $(x, k), A{k), which we defined by equations 
(2.10), (2.11). 

Example Let u{x),—oo < x < oo, he defined by equation (1.25). It is well known that in this 
case 

k-ip ^ ,, k + iptanhp{x - xo) ji.™ , a 

where B{k) is defined by equation (A. 10) with (p replaced by <I>. On the other hand, if we restrict 
u{x) to the half line, it follows that il}{x, k) is still given by (A. 13c) but a{k), b{k), are given by the 
more complicated formulae 

MO,k)+ik^p{0,k) Vx(0,-/c) + i/cV'(0,-A;) 

"^^^ = 2zfe ' ^^^^ = 2rk • 

Thus if u{x) is given by (1.25) it is more convenient to use the completeness relation of the full line. 
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